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We calculate the transverse and the time-time components of the Coulomb-gauge gluon propa- 
gator in SU(3) lattice Yang-Mills theory both on isotropic and anisotropic lattices. The problem 
of scaling violation observed on the isotropic lattice is drastically reduced as the anisotropy in- 
creases; namely, the system approaches the Hamiltonian limit. In the infrared region, the trans- 
verse gluon propagator exhibits a turnover and the temporal gluon propagator shows divergent 
behavior. 
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1. Introduction 

Coulomb gauge provides a very clear picture of color confinement. Coulomb gauge is a physi- 
cal gauge in the sense that the color Gauss' law can be formally solved, and only transverse degrees 
of freedom appears as dynamical degrees of freedom. The striking feature of Coulomb gauge is that 
an instantaneous interaction shows up in the Hamiltonian, which is requisite for color confinement. 
In the Gribov-Zwanziger scenario, the path integral is dominated by the configurations near the 
Gribov horizon where the lowest eigenvalue of the Faddeev-Popov (FP) ghost operator vanishes. 
||l|]. This results in an enhancement of the near-zero modes of the ghost operator, and it has been 
confirmed by the lattice simulations [Q, ||]. Accordingly, the color-Coulomb instantaneous interac- 
tion becomes a confining interaction. The color-Coulomb potential can be obtained by measuring 
the correlator of the temporal link variable at fixed time, and lattice QCD simulations exhibit that 
the color-Coulomb potential rises linearly at large distances and its string tension is larger than the 
string tension of the static Wilson potential [f|, [5[ |6[ 0], which is expected from the Zwanziger's 
inequality [JsJ] . On the other hand, the color-Coulomb potential has been evaluated by inverting the 
FP ghost matrix, and it has been shown that the color-Coulomb string tension almost saturates the 
Wilson string tension [0]. 

The transverse gluon propagator is expected to be suppressed in the infrared (IR) region due 
to the proximity of the Gribov region in the IR direction in the Gribov-Zwanziger scenario flTc|]. 
The instantaneous transverse gluon propagator has been measured by Monte Carlo simulations 



[11, [12|, 13, 14], and recent studies have revealed that it shows scaling violation [13, 14]; namely, 
the gluon propagator calculated at different lattice couplings does not fall on top of a single curve 
after multiplicative renormalization. 



In order to circumvent the problem of scaling violation, the authors of [ ]13[ ] have measured the 
unequal-time gluon propagator 

D«(p,p 4 ) = (A(p,p 4 )A{-p-p 4 )) (1.1) 

and extracted the equal-time propagator D tr (|/?|) by eliminating the p 4 dependence of the unequal- 
time propagator. It has been concluded that D u (\p\) is multiplicatively renormalizable in the Hamil- 
tonian limit and well fitted with the Gribov-type form of the propagator. The method was applied 
only to the transverse gluon propagator and it was not studied if scaling violation can be solved by 
this procedure. 



In pl4J, a new momentum cut is introduced in addition to the cone cut and the cylinder cut, by 
which high momentum data that suffer from discretization errors are excluded from the analysis of 
the instantaneous propagators. It has been shown that this procedure successfully reduces scaling 
violation for the transverse gluon propagator while it fails for the time-time component of the gluon 
propagator. 

The problem of scaling violation of the instantaneous propagator can be seen even at the tree 
level on a finite temporal lattice spacing (we refer to a forthcoming paper for an explicit calcula- 
tion). The reason is that the energy integral does not run from — oo to oo but from —2/a t to 2/a x 
on a finite lattice, and this introduces the spurious \p\ dependence on the free equal-time propaga- 
tor. Therefore, we expect that the instantaneous propagator is multiplicatively renormalizable in 
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the Hamiltonian limit t, = a s /a % — * °°. To make this point clear, we calculate the transverse and 
temporal components of the instantaneous gluon propagator on anisotropic lattices. 

2. Lattice setup and observables 

The lattice configurations are generated by the heat-bath Monte Carlo technique with the stan- 
dard Wilson plaquette action, 

S =J- L 9fcTr(l-I/y + L 5HeTr(l-^ 4 (n)). (2-1) 

? B n,i< j<3 n,i<3 

Here U^ v (n) indicates the plaquette operator, and /3 = 2N c /g 2 is the lattice coupling. On the 
isotropic lattice, the bare anisotropy £g is 1 and the action can be written in a familiar form 

S = f3 £ <fteTr(l-t/ MV (n)). (2.2) 

%b differs from the renormalized anisotropy E, which is defined as the ratio of the spatial lattice 
spacing to the temporal lattice spacing. The ratio of ^ and E, can be determined non-perturbatively 
by matching the spatial and the temporal Wilson loop on anisotropic lattices. We use the relation 
obtained by Klassen for the range 1 < £g < 6 and 5.5 < /? < oo [|l5|]; 

L 1 + (l-i)f^, <2.3> 

where ao = —0.77810, a\ = —0.55055, and 

1.002503gj+ 0.39100^ + 1.47 130^- 0.19231 
T? ^~ Q + 0.26287^| + 1.59008^-0.18224 ' 



We adopt the values of the lattice spacing given in [16] for ^ = 2 and in [17] for ^ = 4, where the 



static quark potential was measured to set the scale. For the isotropic lattice, the scale is set by 



using Necco-Sommer scaling relation Jl 8Q . In our simulations, the first 5000 sweeps are discarded 
for thermalization, and we measured the equal-time gluon propagator for 100 configurations, each 
of which is separated by 100 sweeps. All the lattice parameters are given in Table |. 
In Coulomb gauge the transversality condition 

diAi(x,t)=0 (2.5) 

is imposed on the gauge fields at each time slice, where i runs from 1 to 3. On a lattice, gauge 
configurations satisfying Coulomb gauge condition can be obtained by minimizing the functional 

F u [g] = '£Y i mtTr(l-U?(x,t)), (2.6) 

i=\ x 

defined on each time slice. Here Uf(x,t) = g(x,t)Ui(x,t)g^ (x + i,t) is the gauge-rotated config- 



uration. The functional derivative of Eq. ( |2.6[ ) with respect to g reproduces the Coulomb gauge 
condition in the continuum limit with the linear definition of the gauge field, 



2iga 



(2.7) 

traceless 
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% = a s /a T 


L 3 x L z 


P 
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aj l [GeV] 


a s [fm] 


V[fm 4 ] 
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32 4 


5.70 
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1.160 


0.1702 


5.45 4 


100 




48 4 
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8.17 4 


40 




32 4 


5.80 






1.446 
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4.37 4 
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100 
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4.47 4 


40 
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6.20 
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2.17 4 
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16 3 x 32 
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1.104 
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2.86 4 
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4.29 4 


100 




16 3 x 32 


6.00 
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0.1227 


1.96 4 


100 
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24 3 x 48 












2.94 4 


100 




16 3 x 32 


6.10 


1.718 


1.889 


0.1045 


1.67 4 


100 




24 3 x 48 












2.51 4 


100 




16 3 x 64 


5.75 


3.072 


1.100 


0.1794 


2.87 4 


50 




24 3 x 96 












4.31 4 


50 




32 3 x 128 












5.74 4 


50 




48 3 x 192 












8.61 4 


50 




16 3 x 64 


5.95 


3.159 


1.623 


0.1216 


1.95 4 


50 




24 3 x 96 












2.92 4 


50 


4 


32 3 x 128 












3.89 4 


50 




48 3 x 192 












5.84 4 


50 




16 3 x 64 


6.10 


3.211 


2.030 


0.0972 


1.56 4 


50 




24 3 x 96 












2.33 4 


50 




32 3 x 128 












3.11 4 


50 




48 3 x 192 












4.67 4 


50 



Table 1: Simulation parameters to calculate the equal-time gluon propagator. 



The Coulomb gauge fixing has been done using iterative method with the Fourier acceleration [|TJ], 
and the gauge fixing is stopped if (d,-A,) 2 < 10~ 14 at each time slice. 

We calculate the transverse and time-time component of the equal-time gluon propagator, 



rfv(x-y) = {Al(x)A b v (y))=Df v (x-y), (2.8) 



in the momentum space, 



Df?(p) = 8 ab Uj-^l)D«(\p\) (2.9) 
D? 4 (p) = 8 ab Z *[\^\ (2.10) 



The gauge field is defined as Eq. (2.7). The dressing function Z 44 is constant for \p\ at the tree 



level. In the Gribov-Zwanziger scenario, this is expected to diverge in the IR limit resulting in the 
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confining behavior of the color-Coulomb potential, which is necessary condition for color confine- 
ment in Coulomb gauge QCD. The transverse gluon propagator is expected to be suppressed in the 
IR region due to the proximity of the Gribov region JTo|] - 



3. Simulation results: transverse gluon propagator 
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Figure 1: The equal-time transverse gluon propagator on the isotropic lattice (top left), on the anisotropic 
lattices with % —1 (bottom left), and with % = 4 (top right). The results for the isotropic lattice and the 
anisotropic lattice with | = 4 on large lattice volume are drawn together in one figure for direct comparison 
(bottom right). The cone cut and the cylinder cut are applied and the propagator is renormalized to unity at 
p — 2 [GeV]. 



The top left panel of Fig. [l] shows the transverse gluon propagator on the isotropic lattice at 
various lattice couplings, /3 = 5.7,5.8,6.0,6.2. The cone cut and the cylinder cut are applied [2C] 
and the propagator is renormalized such that D tr {\p\ = 2 [GeV]) = 1. We see that the data points at 
different lattice couplings cross at the renormalization point \p\ = 2 [GeV] and deviate from each 
other both at small and large momenta, as has been observed in 

The simulation results on the anisotropic lattice are drawn in the left bottom {£, = 2) and the 



right top = 4) panel of Fig. |T| On the anisotropic lattice with £, = 2, scaling violation becomes 
moderate compared to the isotropic result, although small deviations among the data points for 
different lattice couplings can be seen. Further increase of ^ leads to a nice scaling behavior and 
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the data points for E, = 4 almost fall on top of one curve, indicating that the equal-time gluon 
propagator is multiplicatively renormalizable in the continuum limit (or in the Hamiltonian limit 
§ — > oo). Accordingly, our results on the anisotropic lattice support our expectation that scaling 
violation observed in the equal-time transverse gluon propagator disappears in the limit £ — > oo. 

In the right bottom panel of Fig. Ill, the instantaneous transverse gluon propagator on the spatial 
lattice extent L = 48 is plotted both for the isotropic lattice and the anisotropic lattice with E, = 4. 
We observe that the propagator has a maximum at p = 0.4 ~ 0.5 [GeV] irrespective of the lattice 
coupling and the anisotropy, and it decreases with the momentum in the IR region. 

4. Simulation results: time-time component of the gluon propagator 




p [GeV] p [GeV] 



Figure 2: The dressing function of the time-time component of the gluon propagator on the isoUopic lattice 
(left), and on the anisotropic lattice with E, = 4 (right). The cone cut and the cylinder cut are applied and the 
dressing function is renormalized to unity at p — 2 [GeV]. 

The dressing function of the time-time component of the gluon propagator is shown in Fig. || 
for the isotropic lattice (left panel) and the anisotropic lattice with E, = 4 (right panel). On the 
isotropic lattice, Z 44 (|/?|) shows scaling violation and the deviation of the two curve is pronounce in 
the IR region. Although the Gribov-Zwanziger scenario predicts that the temporal gluon propagator 
diverges stronger than the simple pole \/\p\ 2 , the numerical result shows that it bends down at small 
momenta for j8 =5.7. 

On the anisotropic lattice, the dressing function shows a much better scaling behavior than 
that on the isotropic lattice. Although the small deviation can be seen both in the IR and ultraviolet 
region, one can expect that the scaling behavior is completely recovered in the Hamiltonian limit. 
Moreover, we find that the IR behavior of Z on the anisotropic lattice is completely different from 
that on the isotropic lattice. For the isotropic case, we see that the dressing function bends down at 
small momenta at /3 = 5.7. By contrast, Z 44 continues to rise with decreasing the momentum even 
for the coarsest lattice data (j8 = 5.75), and Z 44 at available smallest momentum for the anisotropic 
case is about 10 times larger than that for the isotropic case. We note that the spatial lattice spacing 
for (£,j3) = (4,5.75) is larger than that for (£,j3) = (1,5.70). This implies that Z 44 is very sensitive 
to the disretization effects, and taking the Hamiltonian limit is crucial to cure scaling violation for 
the temporal gluon propagator and to explore the IR divergent behavior in Coulomb gauge QCD. 
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5. Summary and conclusion 

We calculate the transverse and time-time components of the equal-time gluon propagator 
both on the isotropic and the anisotropic lattices. We find that scaling violation observed on the 
isotropic lattice is drastically reduced by calculating the propagator on the anisotropic lattices, i.e., 
by getting close to the Hamiltonian limit. In the IR region, the transverse gluon propagator is 
strongly suppressed and shows the turnover at about 500 [MeV]. The time-time gluon propagator 
on the anisotropic lattice is much more enhanced in the IR region compared to that on the isotropic 
lattice. 
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